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ABSTRACT 

Coherent  Anti-Stokes  Raman  Spectroscopy  (CARS)  is  a  new 
type  of  Raman  Spectroscopy.   The  phenomenon  is  associated  with 
the  nonlinear  conversion  of  two  laser  beams  in  a  medium  into  a 
third  collimated  beam  at  the  anti-Stokes  frequency  of  the 
medium.   This  analysis  evaluates  the  performance  capabilities 
of  CARS  under  various  experimental  conditions  for  the  molecules 
Op/N^,H_  and  CO.   Laser  intensity  fluctuations  and  shot  noise 
are  introduced  as  noise  sources  to  predict  signal-to-noise 
ratios  (S/N) .   The  S/N  is  evaluated  as  a  function  of  the  par- 
tial pressure  of  a  gas  to  measure  the  performance  of  CARS.   An 
analysis  of  the  uncertainty  in  measuring  the  rotational  tempera- 
ture is  presented  for  diatomic  gases.   Inverse  bremsstrahlung 
in  plasmas  is  investigated  as  a  possible  process  by  which  the 
parameters  measured  by  CARS  might  be  altered.   An  experiment 
is  discussed  in  which  an  attempt  was  made  to  remove  the  back- 
ground signal  inherent  in  most  CARS  experiments.   Partial 
interferometric  cancellation  of  two  CARS  signals  was  observed 
by  the  proper  positioning  of  two  sample  cells  in  an  otherwise 
standard  CARS  experiment. 
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I.   INTRODUCTION 

Coherent  Anti-Stokes  Raman  Spectroscopy  (CARS)  is  a  method 
of  chemical  analysis  yielding  molecular  spectroscopic  informa- 
tion by  utilizing  nonlinear  scattering  phenomena.   As  a  part 
of  the  growing  field  of  nonlinear  optics,  CARS  has  been  used 
to  measure   nonlinear  coefficients  in  solids,  liquids, 
and  gases  [1-14].    Much  of  today's  efforts  are  focused 
on  the  use  of  CARS  as  a  nonreactive  probe  in  environments  not 
conducive  to  other  techniques  for  spectral  analysis  [14] . 

The  intent  of  this  thesis  is  to  assess  the  usefulness  of 
CARS  as  an  analytical  tool  using  the  example  of  diatomic  gases. 
CARS  is  related  to  the  more  basic  phenomenon  of  Raman  Scatter- 
ing (RS) .   Raman  scattering  involves  a  shift  in  the  frequency 
of  a  small  portion  of  the  photons  incident  on  a  Raman  active 
medium.   RS  has  been  utilized  extensively  to  obtain  molecular 
spectroscopic  information.   CARS  yields  the  same  type  informa- 
tion as  RS  except  that  the  collimated  CARS  signal  is  much 
easier  to  detect. 

CARS  was  feasible  for  chemical  analysis  only  after  the  ad- 
vent of  tunable  lasers  of  moderately  high  power.   The  initial 
worlc  with  CARS  was  performed  by  MaJcer  and  Terhune  [15]  in  1963, 
but  only  recently  has  the  field  begun  to  expand  rapidly,  as 
evidenced  by  the  large  number  of  published  articles.   To  per- 
form CARS  in  gases  requires  lasers  of  much  higher  power  than  in 
liquids  and  solids.   Typically,  pulsed  lasers  on  the  order  of 
megawatts  are  used. 


The  use  of  such  high  power  pulsed  lasers  has  an  inherent 
problem  in  that  the  power  output  from  pulse  to  pulse  is  not 
constant.   This  introduces  a  source  of  noise  that  hinders  the 
observation  of  weak  resonances.   In  addition,  in  a  mixture  of 
gases,  a  CARS  signal  from  other  than  the  gas  of  interest  is 
present  as  a  background  signal.   These  two  factors  represent 
at  present  a  source  of  interference  in  observing  weak  molecular 
resonances. 

As  an  analytical  tool  CARS  can  yield  information  about 
molecular  identity,  concentration  and  rotational  temperature. 
Numerical  estimates  for  the  performance  capabilities  of  CARS 
are  based  on  the  ability  to  calculate  the  third-order  nonlinear 
dielectric  susceptibility,   x    •   This  susceptibility  is  com- 
posed of  a  resonant  term,   y  ,   and  a  nonresonant  term,   Xx-,o  * 

R  NR 

X„   is  calculable  from  the  physical  equations  that  define  it 

but  XvTTj  ^ss  deduced  from  experimental  observations. 

One  source  of  observations  for  x^tt^   relied  on  results  from 

NR 

an  experiment  involving  Stimulated  Raman  Scattering  (SRS) . 
Appendix  A  contains  a  development  of  the  equations  related  to 
SRS  and  shows  the  relationship  of  the  SRS  gain  equation  to 
CARS. 

A  particularly  difficult  environment  to  perform  convention- 
al spectral  analysis  with  Raman  scattering  is  in  an  ionized 
gas  called  a  plasma.   CARS  is  being  used  in  plasmas  and  is 
producing  better  results  than  conventional  Raman  scattering 
techniques.   An  analysis  is  made  to  predict  the  extent  of  in- 
teraction between  a  plasma  and  the  laser  used  for  a  CARS 
experiment  in  the  plasma.   In  particular,  the  inverse 
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bremsstrahlung  process  for   the  transfer  of  energy  is  ana- 
lyzed. 

An  experimental  attempt  was  made  to  improve  the  basic 
CARS  experiment  by  the  elimination  of  the  background  signal. 
A  description  of  the  experiment  and  the  results  are  included 
in  this  study. 


II.   NATURE  OF  THE  PROBLEM 

The  intent  of  this  study  was  to  theoretically  investigate 
the  analytical  capabilities  of  CARS.   The  only  area  investi- 
gated was  concerned  with  diatomic  gases.   It  was  desirable  to 
estimate  how  small  a  concentration  of  gas  could  be  detected 
and  with  what  accuracy  the  rotational  temperature  could  be 
measured.   CARS  has  demonstrated  its  usefulness  as  a  tool  for 
spectral  analysis.   It  offers  orders  of  magnitude  greater 
sensitivity  than  Raman  scattering  for  low  pressure  gases.  [16] 

Typically,  CARS  experiments  are  performed  by  focusing  two 
pulsed  laser  beams  of  different  frequencies  into  a  sample  cell. 
The  beams  cross  at  an  angle  that  corresponds  to  the  momentum 
conservation  requirement.   For  gases  this  angle  is  sufficient- 
ly small  that  the  experiment  may  be  performed  with  the  beams 
colinear.   One  of  the  lasers  is  tunable  so  that  the  difference 
in  the  frequency  of  the  two  may  be  changed.   When  the  differ- 
ence in  the  frequency  becomes  close  to  an  appropriate  natural 
molecular  frequency,  the  CARS  signal  increases  markedly  in 
amplitude. 

The  intensity  of  the  CARS  signal  is  dependent  on  the  pro- 
duct of  the  square  of  the  intensity  of  one  beam  and  the  first 
power  of  the  intensity  of  the  other.   Any  amplitude  instability 
in  the  beams  leads  to  amplitude  instability  in  the  CARS  signal. 
Also,  the  presence  of  an  undesirable  background  signal  masks 
very  small  signals  of  interest. 
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A  brief  sketch  of  several  publications  in  the  field  of 
CARS  has  been  included  in  this  section  of  the  study.   The 
basic  equations  related  to  CARS  have  also  been  included  to 
develop  the  framework  on  which  the  numerical  calculations 
were  based. 

A.  BACKGROUND 

Maker  and  Terhune  [15],  in  1963,  were  the  first  to  perform 
a  CARS  experiment.   Prior  to  the  development  of  good  tunable 
laser  sources  most  of  the  work  with  CARS  was  developed  in  order 
to  measure  nonlinear  coefficients  of  materials. 

Yablonovitch,  Bloembergen  and  Wynne  [17]  were  among  the 
first  to  use  a  tunable  laser  to  generate  a  continuously  tunable 
difference  frequency.   They  investigated  solid  state  materials 
such  as  InSb  and  GaAs .   Levenson,  Flytzanis  and  Bloembergen 

[1]  followed  closely  with  an  experiment  in  diamond  to  measure 

(3) 
the  variations  in   X     as  a  function  of  the  difference  fre- 
quency.  DeMartini  [7]  utilized  CARS  with  H^  gas  in  1972  and 
Taran  [4]  in  19  73  reported  that  he  had  used  CARS  as  a  means  of 
determining  the  spatial  distribution  of  H^  gas  in  a  flame. 
Taran 's  recent  work  has  resulted  in  the  spectra  of  several 
gases,  temperature  measurements  and  spatial  distributions  of 
gases  in  flames  [14] . 

B.  NONLINEAR  EFFECTS  IN  ISOTROPIC  MEDIA 

Due  to  the  nonlinear  dielectric  property  of  materials  two 
laser  beams  may  be  mixed  in  a  medium  to  generate  a  third  co- 
herent beam  such  that 
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oj-  =  2  a),  -  cOp  (1) 

where   o),.   is  the  anti-Stokes  frequency  generated  from  the 
pump  and  Stokes  frequency  (03,   and   002   respectively). 

The  polarization  vector  associated  with  a  material  may  be 
expressed  as  a  power  series: 

where   x     is  the  dielectric  susceptibility  tensor  of  rank 

i+1   and  E   is  the  applied  electric  field.   The  second  order 

(2) 
term   x     is  equal  to  zero  for  isotropic  liquids  and  gases 

because  of  inversion  symmetry.   The  lowest  order  term  that  is 
present  in  a  gas  is  therefore   x    r    the  third-order  nonlinear 
susceptibility.   There  are  several  'physical  processes  that  in- 
volve the  third  order  term,  e.g.,   3co,^aj-.,   2a)-  +  ui^-^u^^, 
2a)-  -  a)_-»'a)_/  etc.   The  last  process,   2a),  -  a)^,  is  the  pro- 
cess referred  to  as  CARS. 

By  assuming  all  fields  are  along  one  axis  the  components 
of  the  tensor,   X/   and  the  electric  field  vector  may  be 
treated  as  scalars.   The  magnitude  of  the  electric  field  may 
be  expressed  as: 

,      i  (k.  z  -  a),  t) 
E^(a))  =  I  [e^e    ^      ^    +  c.c]  (3) 


where   e.   is  the  amplitude,   k.   is  the  propagation  constant 

(equal  to   a).n./c),   z   is  the  length  along  the   z   axis  and 

3 
t   IS  time  in  seconds.   For  the  CARS  process  the  term  E    in 

2 
equation  (2)   includes   3E^   E   [18].   The  other  terms  are  not 

associated  with  CARS. 
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(3) 
The  resulting  third  order  polarization,   P     is  given 


by: 


p(3)_l 
^    8 


-3  (3).  . 

3x   (-0)^0)^^,00^,-002) 


2   * 


i[  (2k^-k2)z- (20)^-0)2)  t] 


+  C.C, 


(4) 


The  factor  of  3  has  been  a  point  of  inconsistency  among  au- 

(3) 
thors  m  the  past.   The  notation  for   y    (-to  ,o),  ,o)  ,o),)   is 
^  '^      a  b   c   d 

consistent  with  that  of  Bloembergen  [19]  for  the  process 


0)    =  0),   +  0)    +  0)  J  . 

abed 


(5) 


Solving  Maxwell's  equations  for  plane  waves  with  polariza- 

(3) 
tion   P    the  anti-Stokes  gain  equation  is: 


de 

3 

dz 


ITTO). 


2   * 


2cn.  ^1  ^2 


^^    £0  (3x    )e 


(3),_i(2k^-k2-k3)z 


(6) 


In  terms  of  intensity,  where   I. 


g—  I  ^.|  ,  after. integrating 


equation  (6)  over  a  length   L,   the  result  is 


^3  = 


A     2 

4lT     0). 


.3C 


3X 


(3) 


2         2,2 
I,    I2L 


^.     /AkL\ 


AkL/z 


(7) 


where   Ak  =  2k--k2-k_.   The  conversion  efficiency  is  defined 

^3 
as  =—  and  if   Ak  =  0   the  conversion  efficiency  becomes  a  , 

2 

maximum.   Because  of  the  small  dispersion  of  light  in  gas, 

CARS  may  be  performed  with   o)  ,  o)   and   o)   colinear.   For  a 

focused  beam  (assumed  to  be  a  cylinder  of  plane  waves  and 

Ak  =  0)  the  conversion  efficiency  becomes: 
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{^)W 


l2...(^]-(t^:A'W" 


2 

P.  .  (8) 


The  third  order  nonlinear  dielectric  susceptibility,  x  r 
is  composed  of  a  frequency  dependent  resonant  part  (Xn)  ^rid  a 
nearly  frequency  independent  part,  the  nonresonant  susceptibili- 

x'^'  =  Xr  ^  X^R.  (9) 

X„   can  be  expressed  [18]  as: 

_  2Nc    da  V (•\r[\ 

^R      ^     4     do.    ~     2       ;  72       ~~     ~7  '        ^^^^ 

3Tia)2        LU)    -  (03,-0)2)   ~  ^Y('^i"'^2 

When   03,-0)^  =  o)  ,   the  resonant  susceptibility  becomes 
1   2    V  f       J 

4 

2Nc    da  .  1  ,,  ,  . 

where  N   is  the  number  density  of  molecules,   (tq)  ^^  ''^^^ 
Raman  cross  section  for  scattering,  A   is  the  difference  in 
probability  of  finding  a  molecule  in  the  ground  and  first 
vibrational  states,  and  y  is  the  line  width,  full  width  at 
half-maximum  intensity. 
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III.   NUMERICAL  ESTIMATION  OF   x 

(3) 
As  defined  m  Section  II. B.,   x     is  composed  of  two 

parts /  the  resonant  and  nonresonant  susceptibility.   The 

resonant  part  is  associated  with  the  various  Raman  active 

frequencies  that  characterize  a  molecule.   The  equation  for 

X„   has  been  derived  [18]   from  considerations  of  a  driven 

harmonic  oscillator  and  quantum  mechanics.   Xtvt-d   ^^  associated 

with  the  electronic  structure  of  the  molecule  and  is  nearly 

frequency  independent. 

Equation  (10)  defines   Xd   i^  terms  of  an  entire  vibra- 

tional  band,  namely  the  ground  state.   To  estimate  Xn      ^^^ 

R 

each  rotational  level  requires  an  analysis  of  each  of  the 
terms  in  equation  (10)  that  are  functions  of  the  rotational 
quantum  number  J.   As  will  be  shown  in  Section  IV,  the  signal- 
to-noise  ratio  for  a  CARS  experiment  is  a  function  of  the  peak 
value  of   Xd*   Being  able  to  predict  signal-to-noise  ratios  is 
important  when  deciding  which  analytical  technique  to  employ 
for  a  certain  situation.   Due  to  the  complexity  of  the  calcu- 
lations a  computer  program  was  utilized  to  predict  the  peak 
value  of  Xo   f°^  several  gases. 

R 

A.   RESONANT  SUSCEPTIBILITY,   Xb 
1.   Cross  Section 

The  Raman  cross  section  for  scattering  is  typically 
quoted  as  the  cross  section  for  an  entire  vibrational  band. 
For  computational  purposes  it  was  necessary  to  apportion  this 
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total  value  among  the  rotational  levels.  The  relative  popula- 
tion of  a  rotational  level  was  used  as  the  apportioning  factor 
such  that 


(da]  ^    (2J+1)  exp[-BJ(J+l)/kT]  (  dg  ] 
dnjj        "^J  Q  \  dfi/, 


(12) 
Total 


where 

Q  =  E   g_  (2J+1)  exp[-BJ (J+l)/kT] , 
J   'J 

J   is  the  rotational  quantum  number,   g    is  the  nuclear  spin 

weighting  factor  and   B   is  the  rotational  constant.   The  values 

used  for  -r^  were  from  Chang  [20] .   His  measurements  of  cross 

section  were  the  results  of  experiments  using  a  ruby  laser  at 

5140A.   The  cross  section  has  a  fourth  power  frequency  depend- 

4 
ence   (w  )  thus  an  accurate  calculation  of  x      requires  that 

the  data  be  corrected  for  this  factor. 

2.  Linewidths 

On  resonance  the  CARS  signal  is  proportional  to  the 
inverse  square  of  the  linewidth.    Two  general  extremes  of 
linewidth  were  used,  the  Doppler  broadened  and  collision 
broadened  cases.   The  Doppler  broadened  lines  are  narrower 
than  the  collision  (or  pressure)  broadened  lines. 

3.  Rotational  Transition  Frequencies 

For  molecules  that  exhibit  Raman  active  vibrational 
and  rotational  frequencies  the  CARS  signal  increases  in  ampli- 
tude as  the  difference  frequency  approaches  one  of  the  molecu- 
lar resonance  frequencies.   The  rotational  frequency  for  a 
vibration- rotation  transition  of  a  diatomic  molecule  is  ex- 
pressed as: 
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0)   ^  =  0)   -2a)x(v+l)-a  J  (J  +  1)  (13) 

V, J    e     e  e  e 


where   co   is  the  haannonic  vibrational  angular  frequency,   x 
is  the  anharmonic  constant,   v   and   J   are  the  quantum  numbers 
of  vibration  and  rotation  respectively,  and  a    is  the  rota- 
tion-vibration interaction  constant. 
4.   Nuclear  Spin  Weighting  Factor 

Homonuclear  diatomic  gas  molecules  have  a  statistical 
correction  factor  for  the  population  of  rotational  levels  due 
to  the  spin  of  each  of  the  nuclei.   Alternating  lines  in  the 
rotational  spectra  have  intensity  ratios   1:(I/I+1)  where   I 
is  the  spin  of  each  nucleus.   The  strong  line  of   D^,0«   and 
N^   are  the  even  numbered  J   levels  while   H-'s  strong  lines 
are  associated  with  odd  numbered  J   levels.   TABLE  1  summar- 
izes the  degeneracy  weighting  factor,   g 


TABLE  1 


NUCLEAR  SPIN  WEIGHTING  FACTORS  (g^) 


Gas 
°.2 

Nuclear 
Spin 

0 

(I) 

Even 

#  J 
1 

Levels 

Odd 

# 

for 

J  Levels 

0 

«2 

1/2 

1/3 

1 

i'2- 

1 

1 

1/2 

Note  that   CO   is  not  a  homonuclear  molecule  and  thus  has  no 
degeneracy  factor;  i.e.,   g.^  =  1   for  all   J. 
5 .   Computer  Program 

To  calculate  x-o      ^^  a  nearly  continuous  function  of 
the  difference  frequency  required  a  computer  program.   A 
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Fortran  program  was  written  (See  page  54)  to  calculate  the 
maximum  amplitude  of  Xn  ^s  a  function  of  the  difference 
frequency  according  to  the  following  equation: 

=2Ncl_  da  A     ^gj^^^-^^)^ ^^r..r     n4) 

^R,  (03,-0)^)  -,+   4  di;^  Q  ^^    2  2  *   ^  ^ 

1  Z      Jho)^         J   0)   T-  -  (0)^-03^)   -iY(w--co_) 


Tables  2  and  3  give  the  results  of  the  computations. 

B.   NONRESONANT  SUSCEPTIBILITY   Xxtt, 

NR 

In  196  7  Rado  [21]  performed  a  series  of  experiments  to 

measure  the  magnitude  of  the  nonresonant  susceptibility,   Xxtd' 

NR 

of  several  gases.   The  impetus  for  the  experiments  was  to  mea- 
sure the  peak  resonant  anti-Stokes-  power  from  a  hydrogen  cell 
and  compare  that  to  the  nonresonant  anti-Stokes  power  generated 
in  a  sample  cell.   The  ratio  of  powers  would  be  equal  to  the 
square  of  the  ratio  of  the  susceptibilities: 


p          /X   (Sample) \^ 
sample     \  NR /     •  H"^^ 

The  value  Rado  chose  to  use  for   Xp   ^as  from  an  earlier 

work  with  phonon  lifetimes  and  SRS  gain  coefficients.   The 

relationship  of   Xp   and  the  SRS  gain  coefficient  is  discussed 

in  Appendix  A.   His  calculated  value  of   x,.tt^   fo^   N^   is 
^  ^NR        2 

LSSxlo"-""^  ^2-  .   Based  on  Table  3,   Xt.   for   N-   at   300  K 
erg  ^  '^P        2 

is  2.4x10    .   Thus,  the  ratio  of  the  peak  to  nonresonant 

erg  ^ 

susceptibilities  is 
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TABLE  2 


3 
Gas        Linewidth  (cm"  )       Temp(°K)     Xp(cm  /erg) 


°2 


"2 


«2 


CO 


Linewidth    (cm 

•1 

.0034 

.0076 

300      8.3  X  10  -"-^ 


1500      2.0  xlO"-"-^ 


.0055  300      3.3  X  10  -^^ 


.012  1500      IJ  X  10  -"-^ 


.036  300      1.0  xio"-*-"^ 


.07  1500      4.3xl0"'-'-^ 


.005  300      3.3  X  10  -^^ 


.07  1500      2.2x10"-^'^ 
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^2 


«2 


CO 


TABLE  3 


Gas       Linewidth  (cm  -*•)       Temp  (°K)     Xp(cm  /erg) 


.16  300     2.6  X  lO"-"-^ 


.07  1500     2.5  X  lO"-"-^ 


.14  300     2.4  X  lO"""-^ 


^2  -17 

^  .06  1500     2.6  X  10 
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.03  300     1.2  X  10  -"-^ 


.08  1500     4.0  X  lO"-*-^ 


.15  300     2.6  X  10  -"-^ 


.011  1500     1.1  X  10  -"-^ 


Xp 

-^  =s  178. 

An  independent  method  to  estimate  the  magnitude  of   X^jp 
was  provided  by  the  spectra  of  N   published  by  Taran  [14] . 
For  the  unresolved  spectra  the  difference  in  frequency  of  the 
maximum  and  minimum  values  of  the  CARS  signal,   5co,   is  given 
by  the  relation 

6a)   =  y|1  +    ^2   •   *  ^^^^ 

Assuming  a  linewidth  of  1.8  cm    and  a  6a)  of  42  cm 
the  calculated  result  is  that 

=s47.  (17) 


By  estimating  the  level  of  the  background  signal  in  the 
N^   spectrum,  another  ratio  of   X-n/X^TD   is  equal  to 


cs  75.  (18) 


Knowing  the  peak  value  of   X^   yields  an  estimate  of   X>td« 
An  earlier  work  by  Taran  [4]  provides  another  means  of 
estimating  the  ratio  of  the  susceptibilities.   The  total 
pressure  of  a   H_,N^   mixture  was  maintained  at  one  atmosphere 
The  anti-Stokes  intensity  generated  was  plotted  as  a  function 
of  the   H-   concentration. 


21 


If   Xp  ^^'^     XvjD   scale  linearly  with  pressure,  then  the 
total  susceptibility  of  the  mixture  may  be  approximated  by 


X^^^  =  XrP  +  tl-PlXNR  <19) 


where   P   is  the  partial  pressure  of  E^-      As  a  ratio  of   x 
to   Y    the  observed  anti-Stokes  power  is  approximated  by 


R 


^3  "Ix'''|^=xL[[^-11  ^'-l]  •        (20) 

When  the  log  of  intensity  is  plotted  against  the  log  of   P 

(See  Figure  1) ,  the  predicted  and  observed  curves  are  in  fair 

4 
agreement  for  a  ratio  of  10  :1.   That  is 


X   2 
^    ^   10^.  (21) 


\r 

Based  on  the  analysis  of  the  H-,N-  plot  and  using  the  value 
of  Xn   fo^   H-   from  Table  3,  the  value  of  x-m-d      ^^^     ^o   ^^ 

estimated  to  be  1.2x10    .   Thus,  another  estimate  of  the 

erg 

ratio  is 


^P    2  4xl0~ 

_£_  =  ^•^■x.-LU     ^  200 

^NR    1.2x10 

The  factor  of  four  between  the  high  and  low  estimates  of  the 
ratio  is  not  surprising  based  on  the  variety  of  sources  used 
to  estimate  them. 
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Figure  1.   Taran's  results  follow  dotted  line.   The  solid 


line  is  the  predicted  curve  for 

3 
two  curves  are  equal  from  10   to 


X 


R  4 

—  =   10    .      The 


^NR 


10      ppm. 


-p 

•H 
CO 

G 
<a 

c 
H  cn 

4J 
W  -H 
Qi  C 
M  D 
O 
•P    >i 

cn  u 
I    <n 

•H   u 

C  -H 
<  i3 

> 

•H 

+J 

IT3 
iH 
0) 
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H^   Concentration  (ppm) 
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IV.   STATISTICAL  CONSIDERATIONS 

The  anti-Stokes  power  resulting  from  a  tightly  focused 
CARS  experiment  is  expressed  as  [18]: 


P3.K 


3X 


(3) 


2    2 
P   P 
1   2 


(22) 


where 


^  ■  (4) 


2    47T  0). 


V    c' 


(23) 


The  resonant  susceptibility  may  be  expressed  in  terms  of  its 
peak  value,   x„^  ^Y- 

IT 


Xr  = 


XpY 


2A(jL)  -  iy 


(24) 


where   Ao)  =  0^,-03^-03  . 
1   2   V 

For  weak  signals  in  which   x  "^^^Xivtd'  ^^^   absolute  value 

squared  of   x     is  given  by 


(3) 


X-dY 


2 


=    Y     + 


NR    2Ao)-iY 


%R  "*" 


2XpYAoj 

2  7    ' 
4(Aoj)  +y 


(25) 


The  magnitudes  of   x     ^t  its  extrema  are   X^piXp/o*   -^^^  ^ 
signal  strength  proportional  to  the  difference  in  the  peak  and 

background  susceptibility 


Xp  2 


NR    %rXp 


(26) 


Thus,  for   Xn  <<  X-KTDf    the  anti-Stokes  power  above  that  due 

P       NR 


to  Xj^R  alone,  is 
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where   rim   is  an  efficiency  coefficient  to  take  into  account 
the  optical  losses,  detector  inefficiencies  and  observed  con- 
version efficiency  factor. 

Typically,  the  largest  source  of  noise  when  performing 
CARS  on  samples  at  one  atmosphere  pressure  is  associated  with 
the  power  fluctuations  of  the  pump  and  Stokes  beams.   If  the 
power  of  the  laser  source,   P  ,   has  a  relative  standard  devia- 

tion  of   £  ,  the  standard  deviation  of  the  power  is  £  P^ .   The 
p  ^         p  L 

fluctuation  of  the  background  signal  power  P^,   due  to   Xvro 
is  thus : 

assuming  each  beam  has  the  same  relative  standard  deviation  z 

For  a  low  intensity  CARS  signal  in  the  absence  of  back- 
ground one  finds  the  major  source  of  noise  to  be  the  shot  noise 
of  the  signal.   For  an  N   photon  signal  the  standard  deviation 
is   /N. 

For  a  pulsed  signal  the  noise  power  due  to  shot  noise  is 

9    P^hv 
P./  =  -^ —  (29) 

N        T 

The  signal-to-noise  ratio  for  the  background  limited  case 
given  by  equations  (27)  and  (28)  becomes 

Xp 

S/N  ^   ^ ,  (30) 

p%R 

whereas  the  signal-to-noise  ratio  for  a  situation  in  which 
there  is  only  a  very  weak  signal  is  given  by  equations  (27) 
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and  (2  9)  to  be 

1       1 
S/N  -  (K  Xp^  ^l^^T^^  (T/hv)^.  (31) 

As  an  application  of  the  signal-to-noise  problem  consider 
the  measurement  of  Xp  from  a  spectral  envelope.  The  uncer- 
tainty of   Xp/  cr     is  expressed  as:   [22] 


■^   =   2.  -^  (32) 


a  2     i  a.2 
P         1 

where  o".  is  the  standard  deviation  in  determining  cr  from 
the  i  pulse  of  a  pulsed  laser.  For  CARS  with  a  background 
limited  situation  and  Xp<<X„p  the  uncertainty  in  the  signal 
amplitude  is  approximately  a  constant ,   cr  ,   associated  with 

the  laser  power  fluctuations. 

9Xp 

Thus  o,    =   G      ( )   where  a        is  basically  the 

1    m   >  m  ^ 

""■^(o).  J,T) 
noise  power  given  by  equation  (28) : 

I,    -r  m\   is  the  intensity  of  the  spectral  envelope  at   co.  , 
(o).  J,T)  J  t'  '^        i' 

due  to  the  quantum  rotational  number  J,  at  temperature   T. 
The  result  of  this  analysis  yields   [16] : 


P    KP  V    X     "^^ 
1     %R 


or 


1      ,2A(jO.  2  /  -DC  N 

^P  =  sTn  ^IfT"^  (25^ 
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In  general  any  standard  deviation  is  dependent  on  the  number 
of  times,   N,   the  measurement  is  taken  such  that 


1 
a   a  — 


TT- 


The  conclusion  is  that  for  a  given  S/N  the  term   (oT^)   is  the 
effective  niimber  of  pulses  required  at  intervals  of  Aoa  to  pro- 
duce an  uncertainty,   a  ,   in  the  measurement  of   Xp* 

The  uncertainty  in  determining  the  rotational  temperature 
is  a  product  of  the  nearly  constant  background  uncertainty, 
a  /   and  the  small  change  of  intensity  with  temperature.   A 

convenient  definition  for   a   is: 

in 

m    max   I 

max 

where   5     is  the  relative  uncertainty  in  the  amplitude  of 

max 
the  line  of  maximum  intensity.   The  uncertainty  in  the  tempera- 
ture then  becomes  [16] : 


^2    ^  2   "J 
^T     '^m 


(!Vt\' 

y  9T    /  • 


^    Z^l  -.^^^    I  .  (37) 


The  intensity  of  a  line  is  proportional  to  the  Boltzmann 
population  distribution.   The  change  of  intensity  with  tempera- 
ture has  been  shown  [16]  to  be : 


^•^J,T  _  ^    [BJ(J  +   1)    -   kT]  ,^P^ 

^T-  -  ^J,T ^^2 •  (^^) 

The  resultant  standard  deviation  in  temperature  is 


27 


T^  [5. 


max 


(39) 


"J,T 


'JJI 


max 


BJ(J  +  1 
kT 


-  1 


It  has  been  estimated  [16]  that  for  an  uncertainty  in  the 
maximum  intensity  equal  to  2%  the  temperature  of   CO   at  1500**K 
may  be  determined  to  within  10 "K. 

This  uncertainty  in  temperature  is  for  an  unrealistic 
situation  in  that  each  rotational  line  is  completely  re- 
solved,  A  more  realistic  situation  is  the  case  of 
a  continuous  envelope  as  given  in  equation  (14).   A  computer 
program  (see  page  59)  was  written  to  calculate  the  uncertainty 
in  a  temperature  measurement,  a    ,      based  on  equation  (37). 


T      ,        31(0).  J,T)  2 
^.  ( k^ ) 


2      2    i 
T     m 


aT 


m    p  p     ^P   p 


1    E 


2      2    i 
T      p 


L  ^P 


.   3(K  X  (i^.  J,T) 

9T 


_   S 
2    i 


X  (w.  J,T) 

±_r 

2 

Xo 


For  the  same  conditions  used  to  calculate  the  previous 
uncertainty  in  temperature  the  results  of  the  envelope  calcu- 
lation yields  an  uncertainty  of  15°K.   This  value  is  of  the 
same  order  of  magnitude  as  the  statistical  prediction  of  10°K, 
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estimated  by  using  equation  (39).   It  is  to  be  expected  that 
the  uncertainty  would  be  higher  for  measurements  on  the 
spectral  envelope.   Verification  of  the  statistical  approach 
means  that  the  much  simpler  method  for  predicting  the  uncer- 
tainty is  a  good  approximation. 
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■  V.  CARS  in;  a  plasma 

Low  pressure  ionized  gases,  or  plasmas,  represent  a  class 
of  materials  that  are  receiving  very  close  investigation.   The 
population  of  the  various  energy  levels  reveals  much  about  the 
plasma.   Determining  the  population  of  levels  is  a  very  diffi- 
cult task  with  conventional  Raman  spectroscopy.   CARS  has  been 
utilized  to  obtain  spectra  of  several  plasmas  at  the  Naval  Re- 
search Laboratory,  Washington,  D.  C.   The  experience  there 
indicates  that  the  technique  is  not  difficult  and  produces 
useful  information. 

Of  importance  when  analyzing  CARS  in  a  plasma  is  how  much 
the  high  power  laser  beams  may_  be  expected  to  perturb  the 
system  under  investigation.   The  major  source  of  interaction 
between  the  plasma  and  laser  beams  appears  to  be  from  the  in- 
verse bremsstrahlung  process  [23]  .   The  free  electrons  in  the 
plasma  will  gain  kinetic  energy  by  this  process. 

If  the  laser  frequency  is  above  the  plasma  frequency  the 
light  will  propagate  through  but  there  will  be  some  attenuation 
or  absorption.   The  plasma  frequency  is  given  by  [24] 


-3 
where   n    is  the  number   density   of  electrons   (cm   )  ,   e 

is  the  charge  of  an  electron  (esu)  and   M   is  the  mass  of  an 

electron  (grams) . 

An  estimate  of  the  change  in  the  electron  kinetic  energy 

may  be  made  by  assuming  there  is  an  absorption  coefficient,   a, 
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for  the  plasma.   The  energy  absorbed  per  unit  volume  of  the 

plasma  is 

I  (1  -  e"°^) AtA   P  Ata 
^ABS  =  ^ S -^A-  <41' 

where   I    is  the  incident  intensity,   P    is  the  incident 
o  .      o 

power.   At   is  the  pulse  duration,   A   is  the  beam  cross  sec- 
tion and   X   is  the  path  length  of  the  medium.   If  all  the 
absorbed  energy  is  transferred  to  the  kinetic  energy  of  the 
electrons  then  the  energy  absorbed  per  electron  is 

aP  At 

e 

As  a  fraction  of  the  initial  energy,   3/2kT  , 


E^„^      aP  At 
ABS        o 


3/2kT     3/2An  kT 


(43) 


where   T    is  the  electron  temperature.   The 'absorption  co- 
efficient,  a,   for  inverse  bremsstrahlung  has  been  reported 
as  [25] 

1     ^2  6 

-   =  ^/^  (375kT)'  yi   3/2  3  fl  -  -P(-|t^1       (^4) 
'       h  c  M  ^  V 

e 

where   n.   is  the  ion  number  density,   z   is  the  charge  of  the 
ions  and  v  is  the  frequency  of  the  incident  light.   Evalua- 
tion of  equation  (43)  with  a   defined  by  equation  (44)  results 
in  a  relative  change  of  temperature  equal  to 

9T     1.3xl0^^n.z^P  At[l  -  exp(-hv/kT)] 

^  =  i 2 (45) 

T  3,  ,,^  ,3/2  •      ^^^' 

e  V  A(kT  )  ^ 

e 
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This  agrees  favorably  with  an  equation  for  the  change  of 
electron  temperature  given  by  [23] 


25    2  T 

AkT     1.4x10   n.z  [1  -  exp(-hv/kT)]/  P^dt 

(46) 


e  _  'i         "^     '     *b    L 


e 

where   (kT  )   is  in  electron  volts  and   P    is  the  power  of 
the  laser  in  watts. 

The  equipartion  time,   t  ,  is  the  average  time  required 
for  the  electrons  to  transfer  their  excess  kinetic  energy  to 
the  ions  and  thus  reach  an  equilibrium  temperature.   The  equa- 
tion for   t    is  given  by  [25] 

25.2  A(T  )^^^ 
t^  ^   r^ (47) 

where  A   is  the  atomic  weight  of  the  ions,   T    is  the  electron 
temperature  in  degrees  Kelvin,   and   ri    is  in   cm   .   If  the 
duration  of  the  laser  pulse  is  shorter  than   t   then  any 
energy  the  electrons  gained  by  inverse  bremsstrahlung  will 
not  be  transferred  to  the  ions  and  there  will  be  very  little 
perturbation  of  the  molecules  under  investigation. 

Consider  the  following  as  a  typical  example: 

4 
plasma:   Deuterium  at  50  Torr,Te  =10   K 

11    -3 
^e  =10    cm   ,  20  cm  long 

-4   2 
laser:    5500A,  5MW,  10   cm  beam  area,  20  nsec 

duration 

The   plasma    frequency,    co    ,    equals    1.8    x   10         sec        and   txi    >>ijd    . 

The  equipartion  time,  t  ,  equals  20  micro  seconds.   The 
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relative  change  in  electron  temperature  given  by  equation  (46) 
is  .0001%.   Little  perturbation  of  the  plasma  is  expected  from 
these  conditions. 
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VI.   BACKGROUND  SUPPRESSION  EXPERIMENT 
As  discussed  earlier,  the  nonresonant  susceptibility,  Xr.iT,' 

NR 

is  responsible  for  the  background  signal.   The  power  amplitude 
fluctuations  of  the  lasers  produce  noise  on  this  signal.   The 
resonant  signal  above  the  background  also  has  noise  on  it.   It 
would  be  much  easier  to  detect  the  resonant  signal  if  there 
were  no  noise  on  the  background  or,  better  yet,  no  background 
at  all.   The  desire  for  reduction  of  the  background  signal  led 
to  an  idea  [26]  for  an  experiment  to  determine  if  it  was  possi- 
ble to  completely  suppress  the  background.   In  addition  to 
yielding  much  better  signal  quality  the  elimination  of  the 
background  would  limit  the  detectibility  of  CARS  to  the  limit 
of  conversion  efficiency. 

Typically,  the  diluent  of  a  sample  provides  the  largest 
contribution  to   X>jp*   "^^^  basic  idea  for  suppression  of  the 
background  signal  was  to  perform  two  CARS  experiments  in  series, 
In  the  first  experiment  the  generated  anti-Stokes  beam,   oo-., 

would  be  due  to  only   y„„   (of  a  diluent)  and  in  the  other 

■^       '^NR 

I 

experiment   the   anti-Stokes   beam,      oa^,    would  be  due   to      Xvrn   "*" 

X„   (the  sample  plus  the  diluent)  .   If  the   oj^   beam,  from  the 

first  cell,  were  then  forced  out  of  phase  by  Ntt   radians 

(where   N   is  an  odd  integer)  relative  to  the  second  beam,  co^, 

and  superimposed  on  the  second  then  the  resultant  signal  should 

be  only  that  due  to   Xo   of  the  sample. 

Molectron  Corporation  of  Sunnyvale,  California,  through 

Dr.  G.  Klauminzer  provided  equipment,  facilities  and  assistance 

for  conducting  the  experiments. 
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Figure  2.   Arrangement  of  Equipment  for  Background  Suppression 
Experiment. 


Figure  2  shows  the  basic  arrangement  found  to  be  most 
convenient.   A  pulsed  nitrogen  laser  was  used  to  drive  two 
tunable  dye  lasers  in  an  arrangement  reported  by  Dr. 
Klauminzer  in  Reference  12. 

Benzene  was  chosen  as  the  first  material  to  work  with 
since  it  was  known  that  the  CARS  signal  was  very  strong  and 
visible.   For  benzene  the  values  for   co,  ,    co^  and  oo-   were 
480.,  504.  and  458.2  nm,  respectively. 

Both  cells  were  filled  with  benzene  in  anticipation  of 

I 

demonstrating  complete  cancellation  of  oo-  with   co_.   It  was 

known  that  a  translation  of  cell  one  along  the  focal  axis 
would  change  the  intensity  of   oo-.   by  changing  the  volume  in 
the  focal  region  but  it  was  also  discovered  that  the  same 
translation  changed  the  phase  of  oo-,   relative  to   oo-.   To 
completely  cancel   oa-.  it  was  known  that  several  conditions 
had  to  be  met.   The  beam  had  to  be  tt  radians  out  of  phase,  the 
intensities  had  to  be  equal  and  the  beams  had  to  be  coinci- 
dent.  A  plot  of  the  resultant  intensity  versus  the  translation 
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of  cell  one  revealed  that  there  was  constructive  and  destruc- 
tive interference  of  the  two  beams.   The  destructive  interfer- 
ence was  not  complete  though.   A  visual  observation  of  the 

I 

0)-  +  CO-   beam  revealed  interference  fringes.   By  steering  the 

I 

00-   beam  off  of  the   co_   beam  the  fringe  pattern  could  be 

changed.   It  was  suspected  that  chromatic  aberration  in  the 
focal  lens  was  responsible.   Calculation  of  the  effect  of 
chromatic  aberration  indicated  that  the  fringe  pattern  should 
be  approximately  that  observed.   The  use  of  an  achromat  for 
the  focal  lens  reduced  the  number  of  fringes  visible  but  it 
was  found  that  a  slight  rotation  of  the  lens  changed  the 
pattern  considerably. 

The  best  measured  ratio  of  minimum  to  maximum  intensity 
as  a  function  of  cell  translation  was  observed  to  be  10:1. 
Since  it  was  the  electric  fields  that  were  interfering,  the 
10:1  ratio  amounts  to  about  a  60%  cancellation  factor.   The 
cancellation  factor  was  arrived  at  by  considering  the  incom- 
plete phase  cancellation  of  the  electric  fields.   The  observed 
intensity  is  proportional  to 

10-,       ie   ^ 
E^  e   -^  +  E2  e   ^    .  (48) 

If  the  amplitudes   E,   and   E^   are  equal,  the  resulting 
intensity  is  proportional  to 

P   10       10     -10       -10 

E^(e   -^  +  e   ^)  (e    "^  +  e    ^)  (49) 


or 

I  a  E^(2  +  2  cos(0^  -  02))-  (50) 
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Complete  constructive  interference,  when   0.  =  Q^,      results 
in  an  intensity  four  times  that  of  one  beam,  and  complete  de- 
structive interference  would  result  in  no  intensity. 

It  is  not  surprising  that  on  the  average  complete  inter- 
ference does  not  take  place.   For  a  small  difference  in   0, 

2 
and   0^   the  observed  magnitude  is  slightly  less  than   4E   or 

slightly  more  than  zero. 

An  observed  9:1  intensity  ratio  is  equivalent  to  a  3.6  to 

2 
.4  ratio  on  the  normalized  scale,   E   =1.   Thus,  the  inter- 
ference produces  a  signal  equal  to  40%  of  what  one  beam  would 
produce,  therefore,  a  cancellation  factor  of  60%. 

Next,  a  three  molar  sodium  benzoate  solution  was  used  in 
the  cells  in  place  of  benzene.   Again  there  was  incomplete 
cancellation.   The  best  observed  ratio  of  intensities  was  5:1. 
It  had  been  predicted  that  sweeping  the  frequency  of   oa^   to 

record  a  spectrum  would  not  alter  the  phase  difference  of   co-. 

I 

and   CO...   The  first  observations  made  while  sweeping   oj^   re- 

vealed  that   co-.   and   w-   went  in  and  out  of  phase  over  a 
period  of  about  10  cm   .   By  changing  the  optical  path  length 
of   0)^   or   cjo,   and   co^   by  use  of  microscope  slides  the 
period  could  be  increased  to  about  60  cm 

The  observed  increase  in  the  period  of  the  oscillation 

verified  that  the  optical  path  length  could  be  controlled  to 

I 

shift   0)^   and   oo.   in  and  out  of  phase.  With  the  appropri- 
ate optical  path  length  the  period  of  the  oscillation  could 
be  adjusted  to  infinity. 

A  spectrum  was  recorded  for  the  three  molar  sodium  benzoate, 
two  cell,  series  cancellation  experiment.   Of  interest  from  the 
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spectra  is  the  interference  of  the  background  with  respect  to 
the  noise  levels  observed.   Figure  3  is  a  portion  of  that 
spectrum  depicting  the  observed  resultant  signal.   When  the 
interference  was  constructive  the  noise  was  much  higher  than 
for  the  destructive  interference. 

As  stated  earlier,  the  reason  for  the  experiment  was  to 
increase  the  sensitivity  of  CARS.  An  analysis  of  the  signal- 
to-noise  ratio  for  a  cancellation  experiment  indicates  that 
the  cancellation  has  to  be  much  more  effective  than  was  ob- 
served in  this  experiment  to  gain  an  order  of  magnitude  in- 
crease in  the  signal-to-noise  ratio. 

2 

Assume   b    is  the  fraction  of  background  intensity  remain- 
ing after  cancellation.   Then,  the  S/N  ratio  is  given  by 

S/N  - 


Xp 


similar  to  the  development  of  equation  (30)  .   An  increase  in 

2 

S/N  of  10  requires   b   =  .01/   or  a  99%  cancellation  factor. 

A  400:1  intensity  ratio  would  correspond  to  such  a  cancella- 
tion factor. 

The  experiment  as  performed  was  very  sensitive  to  the 
alignment  and  quality  of  the  optics  involved.   One  of  the 
major  contributions  to  the  incomplete  cancellation  is  probably 
the  diffraction  limited  ability  of  the  center  lens  to  perfect- 
ly  image  the  initial   oa-.   beam  onto   oo-..   Other  macroscopic 
factors  such  as  turbulence  and  inhomogeneity  of  the  samples 
are  other  obstacles  to  complete  cancellation. 
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VII.   RESULTS  AND  DISCUSSION 

The  basis  for  determining  the  areas  of  application  for  a 
new  analytical  tool,  such  as  CARS,  must  depend  on  recognizing 
the  limits  beyond  which  the  technique  cannot  be  expected  to 
yield  the  desired  information. 

The  results  of  this  analysis  are  presented  as  lower  limits 
of  concentration  that  are  expected  to  be  detectable  by  CARS. 
For  a  mixture  of  gases,  the  detectability  of  one  gas  depends 
on  the  ratio  of   Xd   ^^^   Xt.td'  ^^^   same  as  Taran's  plot  for 
H^   and  N^.   Plots  of  the  ratios  for  other  gases  with  N.   as 
a  diluent  are  presented  in  Figure  4.   The  minimum  detectable 
concentration  is  dependent  on  the  required  signal-to-noise 
ratio  for  a  given  experimental  condition.   The  lower  limit  of 
detection  for  a  pure  gas  at  reduced  pressures  would  be  limited 
to  the  conversion  efficiency  if  not  for  the  photon  statistics 
involved.   From  the  signal-to-noise  analysis,  estimates  of  the 
detectable  lower  pressures  are  summarized  in  Figures  5  and  6. 

Based  on  the  results  of  Section  IV  the  rotational  tempera- 
ture of  a  gas  is  expected  to  be  measurable  to  within  1%.  That 
assumes  the  amplitude  fluctuations  of  the  CARS  signal  observed 
remain  less  than  10%. 

In  a  plasma,  the  inverse  bremsstrahlung  process  is  ex- 
pected to  be  a  relatively  insignificant  perturbation  to  the 
parameters  determined  by  a  CARS  experiment. 

The  experimental  attempt  at  cancellation  of  the  background 
confirmed  the  idea  that  interferometric  cancellation  is  a 
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possible  technique  for  improving  CARS  when  the  background  be- 
comes a  problem.   The  observed  results  were  not  of  the  magni- 
tude hoped  for  but  did  serve  to  verify  the  concept  and  to 
recognize  a  possible  source  of  limiting  performance. 

Turbulence  in  the  samples  and  inhomogeneity  may  be  part 
of  the  reason  that  incomplete  cancellation  occurs.   The  major 
contribution  to  the  lack  of  complete  cancellation  is  probably 
the  result  of  diffraction  limited  optics.   The  first  beam 
cannot  be  refocused  through  a  lens  and  have  precisely  the  same 
geometry  that  it  had  originally. 
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Figure  4.   Log-lo^  plots  of  N»   gas  mixtures  at   1  atmosphere 
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300  K,  Lower  coplot  at  1500  K.   x^td  ^^     ^n    assumed 

TO     3  NK         Z 

to  be  1.2  X  10  -^  ^^2_   at  300  K,  and 

.24  X  10     ^-     at  1500  K. 
erg 


42 


o 

H 

El 


en 

M 

o 
z 
I 
o 

Eh 
I 


2 
U 

H 

cn 


SAMPLE  PRESSURE  (ATMOSPHERES) 


Figure  5.   Signal-to-noise  ratios  for  gas  mixtures  at  30  0  K 
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Figure  6.   Signal-to-noise  ratios  for  Doppler  broadened 

gas  at  SOCK  as  a  function  of  pressure.   Noise 
due  to  photon  statistics.   The  assumed  parameters 
are 


n^  =  .01, 
X  =  5500  A. 


P,  =  5MVI, 
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VIII.   CONCLUSIONS 

This  analysis  presents  a  model  by  which  the  limits  of  the 
information  available  from  CARS  may  be  predicted.   At  STP 

typical  diatomic  gases  should  have  a  signal-to-noise  ratio  of 

-3 
1  at  partial  pressures  of  about  10    atmosphere.   Hydrogen,  be- 
cause of  its  narrow  linewidth  and  population  distribution, 

-5 
should  have  a  signal-to-noise  ratio  of  1  at  about  10    atmos- 
phere.  Pure  gases  at  low  pressure  where  the  linewidth  is 
determined  primarily  by  Doppler  broadening,  may  be  detected 
at  pressures  several  orders  of  magnitude  lower  than  this. 
The  ultimate  detectability  depends  upon  laser  power  and  satura- 
tion effects. 

Simple  relationships  have  been  derived  which  allow  the 
prediction  of  the  uncertainty  in  temperature  of  a  gas  for  an 
assumed  signal-to-noise  ratio. 

In  a  plasma  the  inverse  bremsstrahlung  process  is  found 
to  transfer  insufficient  energy  to  the  plasma  to  significant- 
ly modify  the  parameters  measured. 

The  experimental  attempt  at  background  reduction  by  inter- 
ferometric  cancellation  verified  that  the  principles  were 
correct.   Complete  cancellation,  however,  was  not  observed. 
This  is  thought  to  be  due  to  the  diffraction  limitations  of 
optical  lenses. 

One  area  of  additional  research  and  experimentation  that 
could  further  extend  the  applications  of  CARS  is  that  of  the 
background  suppression.   Another  area  of  improvement  needs  to 
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be  that  of  the  fluctuating  intensities  of  the  high  power 
pulsed  lasers  that  are  available. 
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APPENDIX  A 
STIMULATED  RAMAN  SCATTERING 

A.   DEVELOPMENT 

The  interaction  of  two  electromagnetic  fields  in  a  material 
is  such  that  different  combinations  of  the  two  may  produce 
additional  fields.   The  two  additional  fields  of  interest  here 
are  the  Stokes  and  anti-Stokes  frequencies. 

Raman  scattering  has  been  observed  extensively  and  found 
to  produce  the  Stokes  and  anti-Stokes  frequencies  at  plus  and 
minus  the  vibrational  frequency,   co  ,  of  the  molecule  from 
the  incident  light,   o)^. 

'^Stokes  =  '^l  -  '^v  (^1^ 

^anti-Stokes  =  ^l  "^  ^v  ^^2) 

Stimulated  Raman  scattering  predicts  a  gain  in  the  strength 
of  a  field  as  it  traverses  a  media.   The  media  of  interest  in 
this  study  are  gases.   The  gain  will  be  related  to  the  non- 
linear susceptibility  of  those  gases. 

Beginning  with  Maxwell's  equations,  which  contain  various 
assumptions  implied  by  their  form,  the  plane  wave  equation  is 
derived  from  the  curl  of  the  curl  of  the  electric  field,   E. 
Assuming  a  solution  to  the  wave  equation  of  the  form; 

,      -i  (kz  -  oat) 
E  =  J  (e  e  +  c.c. )  (A3) 

where   £   is  the  amplitude,   k   is  the  propagation  constant 
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and  0)  is  the  frequency.   For  simplicity,  all  fields  have 
been  reduced  to  one  dimension  along  the   z   axis.   From  the 
nonlinear  polarization  term  the  various  possible  frequencies 
are  defined  and  their  field  properties  (gain  or  attenuation) 
are  explained  by  computing  the  solution  to  the  wave  equation. 
For  completeness  both  MKSA  and  CGS  results  are  presented. 

The  development  in  MKSA  units  begins  with  Maxwell's  basic 
and  constituent  equations. 

VxH  =  J  +  D         D=  e^E  +  P 

V  X  E  =  -B  J  =  aE 

V-E  =  0  P=£X2  +  P 


NL 


B  =  0  B  =  y  H 


The   curl   of      H      becomes 


VxH  =   aE+eE   +  P   =    aE+£E   +   P.^^.  . 

O  INJ-i 


Then  the  wave  equation,  defined  from   VX(VxE)   is 


-v^E  =  -  A  (V  .  B)  =  -u^  ^  (V  X  h: 


or 

<T  A    V^E  „    1 


E  +  £:  E  -  ^^^  =  -  -  P^,^.  (A4) 

e     ye     e  NL 


Using  the  assumed  form  of  the  solution,  equation  (A3) ,  and 

2     2 
the  assumptions  that   a  =  0   and   8  e/Sz    is  very  small, 

then  equation  (A4)  becomes 
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2 

0)  ,  ^i(kz  -  oot)  -  2  2     ... 

"T"  ^^®  ^  ^  «  \  -L  '^  c   ,   i(kz  -  tot)  ,      , 

2  +  c.c.)  +  ^  (ee  +  c.c.) 

2n 


kc! 
2 

n 


+  c.c.  =  -  —  P 

J      £   NL 


.  II  ei^'^^  -  '^t)  ^  ^^^^1  ^  _  ^p^^^^  (^5j 


The  usual  definitions: 


2 

= 

1 

^ 

^o^o 

k 

/£/£ 
O 

0)  n/c 

have  also  been  used. 

For  the  two  interactive  fields   E,   and  E^   the  nonlinear 
polarization,   P    is  expressed  as:   [18] 

In  Li 

PnL=3v'"e^%-  (A6) 

2 

Three  frequencies  are  possible  from  the   E,  E^   combination; 

(Dp/  203, -ojp  /  2a)^ +(jjp.   If   Ep   is  the  Stokes  field  with  frequency 
(jjp   (assuming   oj,  <  oj^)  then  a  resultant  field,   E_,   of  fre- 
quency  2oa,-a)p  is  the  anti-Stokes  field.   For  the  anti-Stokes 
calculations   P  j.   is  given  by: 

^    -,  *      -i  [  (2k^ -kp)  z- (2oij--a)p)  t] 

^NL   =    ^V    <ri  =  1^2     '^  12+   c.c.)].        (A7) 


The 


resultant    field  may   be      E_      for   the   case  when 

^    -,  *  -i(kpZ    -    a)„t) 

^NL    =    ^V     <ri^l    ^2'  [^  +C.C.].  (A8) 
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1.   stokes  Gain  Equation 

To  solve  equation  (A5)  for  the  Stokes  condition,  equa- 
tion (A8),  there  are  two  possible  alternatives  that  lead  to  a 
gain  equation.   One  way  is  to  assume  that   8£^/9z  is  very 
small  and  k-   is  complex  such  that 

k2  =  k'  +  ik". 

The  imaginary  part  of  k^/  k" ,  is  then  evaluated  as: 

^o  3x"|E  °|2 
k"  =  -      ^ 


2cn 

where   x"  is  the  imaginary  part  of   x   and  E,    is  the 
magnitude  of  E,   at  z  =  0 . 

The  observed  intensity  of  E^      is  proportional  to  the 

I   1 2 
time  average  of   jE^I    or 

.  — 12    1   * 
where 


U  TT  *    1/    ik'2  -k"z  ^    *  -k'z  -k"z.  ,      ^    , 
2^2^2   =  8^^2^    ^      ^    ^2    ^  ^  )(c.c.) 

-2k"z 

Thus,   |E^|   is  proportional  to   e      .   x"  is  a  negative 

-2k"  z 
quantity  [27]  and  therefore   e       is  an  exponentially 

growing  function  with  distance. 

The  other  possible  way  to  solve  equation  (A5)  is  to 

let 

2   2 

,2    ^2  ^2 
k. 


•2        2 


and  then  solve  for  de^/^z.      The  result  is  that 
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8^2       ^ 


=    il^      3(x'  ^  ix")  |E,°|    ie^l 


and  again  for   x"   negative  there  is  gain  for  increasing   z. 

Both  approaches  result  in  an  intensity  gain  coeffici- 
ent, g,      such  that 

1=1  °e^^ 

where  1^        is  the  intensity  of   I-   at  z  =  0 .   The  result  is 
that 

0)-  2 

g  =  7p^   |3x"l  |E3_°  I  ,  (A9) 

For  cgs   units   the   gain   is 


4-nu}  2 

g  =  -^      |3x"l     \E^    \  .  (AlO) 


2.   Anti-Stokes  Gain  Equation 

Solving  equation  (A5)  for  the  anti-Stokes  condition 
(A7)  results  in  the  gain  equation  for  the  anti-Stokes  field. 
The  only  way  to  solve  equation  (A5)  is  to  assume 

2   2 

k.  =      -^ 


•3       2 
c 

The  result  is  that 


^  -    ^   3(X'  +  ix")  h.^.^o*   e"-^^^  (All) 


dz        T]  ^c         '^  ^    '    8  112 


where   Ak  =  2k^  -  k^  -  k-. 

Again  if   x"  is  negative,  there  is  gain.   In  cgs  units 
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=  47T(^)^^,,  (A12) 


az (cgs)      '  az  'MKSA. 

B.   CORRELATION  OF  VARIOUS  ARTICLES 

The  only  published  values  for  nonresonant  susceptibilities 
of  gases  are  those  given  by  Rado  [21] .   It  was  desirable  to 
check  for  consistency  in  the  equations  used  by  Rado  and  the 
ones  used  to  calculate  resonant  susceptibilities. 

The  basis  for  Rado's  experimental  values  of   Xvjn  ^^^  "^^^ 

(3) 
resonant  value  of   x     fo^   ^9   from  other  experimental 

data  [23] .   The  gain  coefficient  for  Stokes  Stimulated  Raman 
Scattering  provided  Rado  with  the  needed  value  of   x-d   ^^^  ^o 
at  pressures  greater  than  10  atmospheres.   The  gain  coeffici- 
ent was  given  in  [28]  in  cgs  units  as 

k[87r^L  A  Na^]  x       -, 

■^   =  ^r X  10   (cm/W)  (A13) 

I,      £  en  n  V   /  / 

where   L   is  the  local  field  correction  factor,   N   is  the 
number  density  of  atoms,  A   is  the  population  ratio,   x    is 
the  phonon  lifetime,  a   is  the  polarization  coefficient  and 
I^   is  the  incident  laser  intensity. 

An  independent  article  by  Maier  [29]  used  the  following 
definition  for  the  gain  coefficient 

2 

g    _      BttN    cog 

■■■l       iix)^    c^   Av 

where  Av  is  a  linewidth.   Rado's  gain  is  equivalent  to  Maier 's 
if 

A  =  1. 
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L  =  1. 

2 

e  =  n 

and 

T    =  1/2TTC  AV. 

r 

The  experimentally  observed  gain  for   H^   above  10  atmos- 

-9 
pheres  was  stated  as  1.45  x  10    cm/W   for  an  incident  light 

source  of  975  nm.   If  the  gain  coefficient  developed  in  equa- 
tion (AlO)  is  divided  through  by  intensity  then 


(A15) 


^  ^  3271^  oa^  |3x"| 

I^        2  2 

1       c  ri 

Solving  for    | 3x" |   the  result  is 

gcA 
|3x"i  =  1  (A16) 

647T 

2 

assuming   n   =1.   The  value  for   3x"   is  the  same  as  that 

used  by  Rado  to  represent   x"« 

The  gain  coefficients  given  in  two  textbooks  more  nearly 
resemble  the  gain  coefficient  in  equation  (A9) .  Pantell  and 
Putoff  [30]  define   g   in  MKS  units  as 

47T  X"  I. 
g  =  ^ i—  (^17) 

^2^2^1^o^ 

whereas  equation  (9)  may  be  rearranged  as 


47T  |3x"|  1-^ 
^  ^      ^2^2^o^ 


(A18) 


differing  from  Pantell  by  a  factor  of   n-j_-   If   n-j_  =  1/   they 
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become   the  same.   Yariv  [27]  gives  a  voltage  gain  coefficient 
as 

g'  =  — ^   |E|^  |X"I  .  (A19) 

The  intensity  gain,  g,    would  be   2g' .   Thus,  the  intensity 
gain  coefficient  is 

2g'  =  g  =  -%  |E|^  |x"|  =  -%  |X"I  |E|^         (A20) 
^2  2 
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